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SUMMARY 


An investigation is made By means of the linearized theory of the 
relations Between the supersonic stability derivatives of a wing and 
its reverse obtained by interchanging the leading and trailing edges . 

The analysis includes wings that have a supersonic leading edge, a 
supersonic trailing edge, and either supersonic side edges or nonre- 
entrant subsonic side edges which are noninteracting, that is, neither 
edge lies in the zone of action of the other. The investigation is 
made for thin flat -plate wings and includes steady vertical and longi- 
tudinal motions, steady rolling and pitching, and motions in which the 
wing is accelerating uniformly in the vertical direction. 

The results of the investigation are obtained with respect to body 
axes and indicate that the rates of change of lift and drag with angle 
of attack and forward speed, the damping in roll, and the rate of change 
of lift with a uniform variation of angle of attack with time were all 
equal for the wing and its reverse . Also, with respect to axes of 
rotation that are fixed with respect to the wing configuration, the 
damping in pitch was equal for the wing and its reverse; whereas, the 
rate of change of pitching moment with angle of attack for the wing was 
equal to the rate of change of lift with pitching velocity multiplied 
by the stream velocity for. the reverse of the wing. 


INTRODUCTION 


The theorem of flow reversibility as demonstrated by Von Kerman 
(reference l) states that, within the validity of the linearized theory, 
the wave drag of an uncambered wing at zero lift is the same as that of 
its reverse obtained by interchanging the leading and trailing edges . 
The theorem of flow reversibility was later extended by Hayes (refer- 
ence 2) to apply to the rate of change of lift with angle of attack for 
wing configurations in which all the edges are supersonic. (The term 
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"supersonic edge" as used herein refers to an edge for -which the normal 
component of the stream velocity is supersonic; a "subsonic edge" is 
defined in an analogous manner.) For this type of wing configuration. 
Flax (reference 3) obtained relations between a wing and its reverse 
for the rate of change of lift and pitching moment with angle of 
attack, for the rate of change of lift and pitching moment with pitching, 
and for the damping in roll. A treatment of the flow-reversibility 
theorem for the rate of change of lift with angle of attack has been 
reported in reference k for wings which have a supersonic leading edge 
and supersonic trailing edge and nonre-entrant subsonic edges which are 
noninteracting. A nonre-entrant edge is one for which no part of the 
trailing vortex sheet impinges on another part of the airfoil. Non- 
interacting subsonic edges are such that neither edge lies in the zone 
of action of the other. 

In the present paper, simple relations between the stability 
derivatives of a wing and its reverse in supersonic flow are derived by 
means of the linearized theory. The analysis includes wings that have 
a supersonic leading edge, supersonic trailing edge, and either super- 
sonic side edges or nonre-entrant subsonic side edges that are 
noninteracting. 

The analysis is made for thin flat-plate wings for steady vertical 
and longitudinal motions, for steady rolling and pitching, and for 
motions in which the wing is accelerating uniformly in the vertical 
direction. The analysis for the wings with side edges that act as a 
subsonic trailing edge (subsonic portions of raked-in wing tips) is 
based on the application of the Kutta-Joukowski condition. 


SYMBOLS 


rectangular coordinates 

auxiliary variables used to replace x and. y, 
respectively / - . 

* -v 

Mach coordinates used^in appendix B (see fig. 6) 

source point in Mach coordinate system 

incremental velocities of wing along x- and z-axes, 
respectively; u is positive in flight direction; 
w is positive downward 


*,y,z 

M 

4 

u,v 

u l> v l 

u,w 
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Uo,Vo 


Auq ' 


8u 0 ',8L' 


V 


/ 

disturbance velocity of fluid along x- and z-axes, 
respectively; Uq is positive in free -stream 
direction; w Q is positive upward 

contribution to Uq for steady motions in reverse flow 
of vorticity shed in wake of portion of subsonic tip 
which is raked inward (based on application of Kutta- 
Joukowskl condition) 

partial contribution of sources in an elemental area to 
velocity Uq ' at a point and to elementary lift over 
an area, respectively 

free -stream velocity 


a speed of sound; also, used to refer to a line and, with 

prime, to refer to a line for wing in reverse flow 
(see figs. 2 to 4) 

M stream Mach number (v/a) 

B =\/m 2 - 1 

R = \/( x - £) 2 - B 2 (y - tj) 2 




1 

?x = 


Sx 


t 


r l> T 2 


*t« 




p 

AP 


Mach angle 

disturbance -velocity potential on upper surface of airfoil 


time following disturbance 


time-lag parameters defined by equations (31) ^nd ( 32 ), 
respectively _ _ 1 AKi 


E ^c^o G uoht\ 1R ckaft 

CBA D Ws,T^ 


mass density 


pressure difference between lower and upper surfaces of 
airfoil, positive in direction of lift 
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A(AP),A(AP) ' contribution to AP in uniform accelerated vertical 

motion in direct flow and in reverse flow, respectively, 
of vorticity shed in wake of portion of subsonic tip 
which is raked inward (based on application of Kutta- 
Joukowski condition) 

C,D,E,F,G regions of integration for obtaining contribution to 0, 

(jfx, or Uq (see fig. 2); also used as subscripts for 
limits of integration in appendix A 

A 

c 

s 


combined regions C, D, E 
points in a region (see fig. 2) 

function representing equation for leading edge (ji t.r 
mean chord 
wing area 


b 

7 


i 

a 


wing span 

local inclination of wing tip to x-axis at point of 

intersection of line a or a 1 (fig. 3 or 4) with wing 
tipj positive when tip slopes to right 

point of intersection of line a with side edge in figure 4 

wing angle of- attack in steady flight ( w/v ) 


a 

•p,q. 

D,L 

®wo 

M,I* 


"m 


derivative of a with respect to time 
angular velocities about x- and y-axes, respectively 
drag and lift forces, respectively 
wave drag at zero lift 

pitching and rolling momenta, respectively 


drag, wave drag at zero lift, and lift coefficients, 


respectively 


D ? Dwo t I*' 

fr 2 s 


pitching-moment coefficient f ■ 


M 
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/ Lp \ 

rolling-moment coefficient / j 


♦ 

additional potential associated with 
integrals in equation (C3) 

A(AP) or each of 

T 

additional potential associated with 
equation (C4) 

integral in 

K 

constant used in equation (C3) 


F 

constant used in equation (C4) 


E 

constant used in equation (C6) 


x ' 

auxiliary variable of integration which replaces x 

Subscripts : 



1,2,3 

contribution of component terms of equation ( 3 ) to 5u 0 
and 6L 1 

LE 

refers to wing leading edge 



Primes used with disturbance velocities, forces, moments, angular 
rotations, and stability derivatives refer to wing in reverse flow. 

Whenever a, q, u, p, and a are used as subscripts, a nondimen- 
sional derivative is indicated and this derivative is the slope through 

zero. For example, 


dC m 
, — \ 

> °Lu = 

dC L 


t 

O 

t 



u->0 


C V= 




\2 V 


and Cl^ = 


p->0 




<§ 


-J tt >0 y "t - 


y 
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ANALYSIS 


The analysis is based on the usual assumptions for the lin earized 
theory of supersonic flow, that is, isentropic and nonviscous flow, 
small disturbance velocities relative to the free -stream velocity, and 
a constant velocity of sound. The wing is assumed to be thin and flat 
and to lie in the xy -plane (z = 0). 


STEADY MOTIONS 

Derivation of Basic Equations 


If th 9 foregoing assumptions are made, the velocity potential (ft 
at any point (x,y) in the plane of the wing for steady motion may be 
expressed in the following form (see, for example, reference 5): 

j((x,y) . - ijjy °(^ )aS d 1 (; 


where B = v(x - |) - B^(y - q)^ and where the integration extends 

over the entire region in the plane of the wing that can influence the 
point (x,y). The quantity w Q U,q) is the fluid perturbation velocity 
in the vertical direction at the point (|,q,0). Equation (l) expresses 
the condition that a solution for (ft that satisfies the basic linear- 
ized partial-differential equation of the flow and the required 
boundary conditions may be obtained by superposing elementary source 
solutions of strength proportional to w 0 over the region of 
integration. 

The present analysis is limited to plan forms that have a super- 
sonic leading edge, a supersonic trailing edge, and either supersonic 
side edges or nonre-entrant subsonic side edges that are noninteracting. 
A general plan form that fulfills these conditions is shown in figure 1. 
The typical region that can influence a point (x,y) is contained in 
the xy -plane within the Mach forecone from the point. Figure 2 illus- 
trates the typical Mach forecones or regions of influence for two 
points P and Q for a wing with a supersonic leading edge, a super- 
sonic trailing edge, and a streamwise tip. The Mach forecone may 
intersect only the leading edge, as for point P in figure 2, or it 
may intersect the side edge as well, as for point Q in the same 
figure. The treatment of the second situation is the more complicated 
and includes the first situation as a special case . The second 
situation is shown in detail in figure 2. Thus the component regions 
that influence the point Q, are C, D, E, F, and G. It has been 
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shown in reference 5 that, for a thin flat-plat© ving in certain cases, 
the influence of the external field F on the point (x,y) is identi- 
cally canceled by the effect of region G in the wing. This cancella- 
tion is valid under the following conditions: 

1. The vorticity in the wake of the wing cannot influence any part 
of the wing. 

2. The regions external to the wing that are influenced by subsonic 
edges must not interact with each other in such a way as to influence 
any part of the wing (noninteracting subsonic edges) . 

Restriction 1 limits the validity of the cancellation to portions 
of a subsonic side edge which are either in the streamwise direction or 
raked out'. Those portions of a subsonic side edge which are raked in, 
however, will produce a wake in which the vorticity will influence 
portions of the wing. (See fig. 1 or 3(b).) 

If the point (x,y) (fig. 2) lies inboard and out of the tip region, 
the exterior field F disappears and the regions of influence for the 
point are bounded by the wing boundary and the Mach lines from its 
forecone. The component regions that influence the point thus degener- 
ate into two regions which correspond to C and E (fig. 2). These 
component regions C and E also correspond to those that inf luence 
any point in a wing having supersonic side edges that are raked in or 
out. Thus, the supersonic raked tip may be treated simply as a con- 
tinuation of the leading-edge wing boundary. 

Because of this degeneration of the component regions C, D, 
and E, the subsequent analysis refers specifically to points on the 
wing that are influenced by subsonic tips. The results obtained for 
this case can' then be made to apply also to the region between the tip 
Mach cones for wings with subsonic tips and to the entire wing for the 
case of supersonic tips of arbitrary rake by simply inserting the 
condition that the component region D vanishes for these cases . 

For simplicity, the analysis that follows immediately refers to 
figure 2 in which the tip is shown to be straight and in the direction 
of the stream. The general analysis, however, applies to all plan - 
forms in which the cancellation of the region external to the wing is 
valid. 

The pressure at the point (x,y) is proportional to the 
quantity $x( x >y) that may be obtained by differentiating equation (l) 
with respect to x. The differentiation of equation (l) with respect 
to x under the integral sign must be performed with special attention 
to the variable limits for £ and q. (See fig. 2.) If this 
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differentiation is performed, the contribution of each component region 
for ^(x,y) (regions C, D, and E in fig. 2) results in five 
integrals. One of the integrals for each region is. a surface integral; 
vhereas the remaining four are integrals that are to he evaluated along 
the boundary of the region. These five integrals for 0 x (x,y) are- 
given in appendix A. 

As a result of an analysis that is given in appendix A, the fore- 
going integrals for the component regions C, D, and E, which are 
hereinafter denoted by A, may be written in the following form: 

UoU,y) = 0xU>y) 

,tQ(x - l)d£ dt) _ ^W\r wo(l,n)<U _ vq U, y) 

B 3 *Vdx/Ja ■ E B 

( 2 ) 

The notation I designates the "finite part of the infinite integral, " 
a concept introduced by Hadamard (reference 6) . For the cases con- 
sidered in this analysis in which w Q is constant, is any function 
of £, or* varies linearly with T), if the integration is performed 
first with respect to q, it has been verified that the q limits for 
the surface integral along boundaries of the Mach forecone from the 
point must be omitted. More generally, however, the surface integral 
in equation (2) for an arbitrary variation of w D -with q can be 
evaluated as the finite part of the integral by the method introduced 
by Hadamard (reference 6) and elaborated in reference 7 • 



Elementary Lifts for Wing in Direct Flow 
and in Beverse Flow 

Equation (2), when modified for the wake effect contributed by 
portions of a subsonic side edge that are raked inward is the basic 
equation for deriving relationships between the stability derivatives 
of a wing and its reverse for steady motion. In the derivation, the 
concept of a forezone of influence and an afterzone of influence of a 
point is used extensively. Thus, consider the wing shown in figure 3(a)' 
in direct flow. As noted previously, only the sources in region A can 
influence the velocity Uq at P. Thus A is called the forezone of 
point P . Now consider the forezones of type A for various points P 1 
below P, which are located so that the forezones from these points 
include the point P. The point P’ is then found to be limited to the 
area A’ in figure 3(a). Thus the area A' defines the region of 
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disturbance produced by a source element at P. The region A' there- 
fore is called the afterzone of influence of point P, or the afterzone 
of point P. 

A simple relationship is subsequently shown to exist between the 
integrated effect of the sources in region A on the elementary area 
at P in direct flow and the integrated effect of the sources in the 
elementary area at P on the region A in reverse flow. 

In figure 3(a) for direct flow, let the source points be tenta- 
tively designated as and field points as (x,y). If P is con- 

sidered as a field point, then the velocity Uq at P or (x,y) is 
given by equation (2) . The following relationship can be obtained from 
figures 2 and 3(a): 

d^2D = _ l /l - B tan 7 ] 
c3x B\1 + B tan 7) 

where 7 represents the local inclination of the wing tip to the x-axis 
at the point of intersection of line a with the tip. Equation (2) 
for direct flow then becomes 


Uo(z,y) = y) 



Suppose now that the direction of the stream flow is reversed. For 
this reverse flow, it is convenient for the purpose of later develop- 
ments to reverse the notation for source points and field points. Thus, 

Direct flow Reverse flow 

Source point: (|,t|) Source point: (x,y) 

Field point: (x,y) Field point: (|, re- 
consider the wing in figure 3(a) vhen the direction of flow is 
reversed and let P be a field point. Then, relative to the reverse- 
flow direction, the line a ' intersects the side edge at a point at 
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which the local slope of the side edge is inward :(raked-in tip portion) . 
This portion of the wing tip acts as a subsonic trailing edge and the 
vorticity shed in its wake influences the wing. (See fig. 3(b) .) The 
contribution of this vorticity or wake effect to the velocity Uq ' for 
a point on the wing in reverse flow is designated as Auq ' . This con- 
tribution Auq' is determined on the basis that the Kutta - J oukowBki 
condition is satisfied along the subsonic side edge. The contribution 
of the part other than the wake effect to the velocity Uq 1 at point P 
or ( !>t|) is obtained from equation (3) in three steps: First, the 

signs of I and x (also 7) are reversed (but not those of the 
element dx) because of the reversal of the direction of the x-axis 
relative to the flow direction. Next, the coordinates |,t) and x,y 
are interchanged in each system because of the changed notation for 
source and field" points in reverse flow. Finally, the area for the 
double integration is changed to A.' and the path of the line integral 
is specified to be along a ' . Thus, the velocity Uq ' at point P 
or (1,1)) in reverse flow is 


ff v °' (i » y)(i - e)dx< fr 

0 it JJ a* r 3 


+ 1 fi ± B tan z ) r v o *(x,y)dx _ w 0 '(g,T)) . , 

*B\L - B tan yj J a * R B ^ 


(k) 


By means of an extension for generalizing the results obtained in 
reference 8, it may be shown that in the case of flat plates in which 
the side edges are raked in and are subsonic, if the Kutta- Joukowski 
condition is applied, the vorticity in the wake is of just sufficient 
strength to cancel the contribution to uo of a line integral that 
yields an infinite pressure along the subsonic tip. This line integral, 
which is canceled by the wake effect, is evaluated on the wing along 
the line a*. Thus, it is found that, with regard to the reverse 
flow shown in figure 3(a), the effective forezone of influence for 
the point P or ( |,t]) is limited to the area A f when the Kutta- 
Joukowski condition is applied, 

Again, by an extension of the results in reference 8, the contri- 
bution to u 0 * of the vorticity in the wake or Au 0 ’ is 


A . 2 tan 7 r w 0 'U,y)dx 

^“° = ' «(1 - B tan 7) Ja' 


R 


(5) 
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Substituting the expression for Auq ' from equation (5) into equa- 
tion (4) and simplifying yields 


Uo' ( M)=i[£ 


wp '(x,y)(x - g)dx dy + _1_ T vo'(x,y)dx _ v 0 '(|,r|) 
A' R3 + «BJa’ R " B 


(6) 


In figure 4(a) is sho-wn a portion of a wing in direct flow and in 
figure 4(b), in reverse flow. The wing is intended to meet the general 
conditions for the plan form illustrated in figure 1. In figure 4(a) 
for direct flow, the shaded region A is the forezone of the 
point (x,y) in the element dx dy. In figure 4(b), for reverse flow, 
the shaded region A is now the afterzone of the same point. Thus, in 
direct flow, the sources in region A induce lift on the element dx dy; 
whereas, in reverse flow, the sources in the element dx dy induce lift 
over region A. For the -cases treated in this paper, a simple relation- 
ship is subsequently shown to exist between these two elementary lifts. 

Elementary lift in direct flow. - For convenience in the analysis, 
let the point (x,y) in figure 4(a) in direct flow be such that the 
line a intersects the side edge at a point i at which the local 
slope of the side edge is outward (raked-out tip portion) . This choice 
does not affect the generality of the subsequent analysis . The lift on 
the element dx dy in figure 4(a) can be evaluated at once. The 
velocity Uq at dx dy is given by equation (3) and the lift 8L 
according to Bernoulli's law is equal to 2pYuo dx dy. Thus, 


8L(dx dy) = 2pV dx dy 


JL 


v 0 (l ,T|)(x - S )dt dT| 

53 


1 Pi - B tan 7 \ P v o( S, 
«B\1 + B tan 7/ J a R 


T))<i£ _ vo(x,y) 

B 


(7) 


Elementary lift in reverse flow.- In the reverse flow shown in 

figure 4(b), the point i is now located relative to this reverse-flow 
direction at a point at which the local slope of the side edge is inward 
(raked-in tip portion). For this case, therefore, equation (6) is used 
to determine the lift on region A due to the sources in dx dy. This 
evaluation is made by subdividing equation (6) into its elementary com- 
ponents, the contributions of which are then resummed in a different 
fashion to give the total lift on the region. 
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Consider first the contribution of the surface integral in equa- 
tion (6) over region A' in figure 3(a) for reverse flow, where 
region A' is the forezone of the point P or (l,q). The partial 
contribution of the sources in is dy (any element in region A') to 
the velocity Uq ' at the point (|,q) is obtained from equation (6) by 
dropping the integral signs. Thus, 


8 lUo'(S,tl) = - 

It 


y 0 '(x,y)(x - l)dx dy 
R 3 


( 8 ) 


In figure 4(b) for reverse flow, therefore, each element d£ dq in 
region A experiences an incremental lift induced by the sources in the 
element dx dy of amount 2pV d £ dq times the quantity given' in equa- 
tion (8). The corresponding incremental lift for the whole of region A 
is then given by the integral 


1 = dx dy 


rr w o’(x,y)(x - £)d£ dq 

JJ A r3 


(9) 


The lift contribution that results from the line integral along a 1 , 
which appears in equation (6), requires special consideration. The 
details are given in appendix 3. The result is 


62L' = 


gpg f 1 - ., B ..t_ a fr -..?W dy 
Bit Vl + B tan 7) 


w 0 ’(^y)dl 


(10) 


This incremental lift is to be regarded as localized in a strip of 
infinitesimal width along the line a. The third part of the contri- 
bution of the sources in dx dy to the incremental lift on region A 
(fig. 4(b)) is obtained from the last term in equation (6). This 

term - v ° expresses the contribution at a point only if the 

field point is in the immediate vicinity of the source point, that is, 
for (x,y) — 5> (l,q). Thus, the point (x,y) in region A experiences 
an Incremental contribution 1 because - of the sources in dx dy of 

amount 


8 3 u 0 '(x,y) = 


w Q '(x,y) 


B 
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The incremental lift due to this contribution is 

& 3 L* = - 2p7v o , (^y) dx ** (11) 


The total lift induced on the region A by the sources in the 
element dx dy in reverse flow is given by the sum of 61 L &gL ', 
and 63!', which are evaluated in equations (9) > (10), and (ll) , respec- 
tively. The total lift is 


BL ' = 2pY dx dy 


1 


v 0 'U,y)U - !)d! drj 
R3 


1 / 1 - B tan p 
*B \1 + B tan 7/ J 


w 0 '(x,y)d£ w 0 '(x,y) 


B 


( 12 ) 


Reversibility Relations for Stability Derivatives 

. 1 

In the subsequent analysis all the stability derivatives are 
expressed with respect to body axes that are fixed with respect to the 
wing configuration. 

In the analysis for the stability derivatives that involve moments 
and angular velocities the following remarks concerning notation are 
pertinent. In conventional notation, the positive direction for moments 
and angular velocities is defined with respect to the direction of 
flight. The flight direction is, in effect, reversed relative to the 
coordinate axes (fig. 4) for the wing in reverse flow. Thus, the sign 
for all moments and angular velocities is also reversed relative to 
these coordinate axes that are fixed in the wing. The positive direc- 
tion for the several moments and angular velocities is indicated in 
figure 5* The analysis has dealt for convenience with the relative 
behavior of a given wing in direct and in reverse flow- However, a 
given wing in reverse flow and the reverse of a given wing in direct 
flow are exactly equivalent and their stability derivatives are 
identical. 

Derivative Ci^.- Reconsider the significance of 6L and BL*. 

With reference to figure 4(a), BL expresses the lift induced on the 
element dx dy by all the sources in the forezone A in direct flow. 
With reference to figure 4(b), BL' expresses the total lift induced on 
the afterzone A by the sources in dx dy in reverse flow. The 
area A is geometrically the same in both direct and reverse flow. The 
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expressions for 8L and 5L' are now recapitulated for comparison: 


6 L = 2pV dx dy 


1_ 

*B 


1 - B tan 7 
1 + B tan 7 


VqU,t|)(x - t)d-l 
B3 

P v 0 ( l,t|)d| w 0 (x,y)j 

a R B J 


and 


8L' = 2pV dx dy 


1 

jt 


w’(x,y)(x - l)di dT) 


l/ l-Btany ) f Vo'(^y)d| _ w 0 '(x,y) 

+ rtB \1 + B tan 7/ J R . B 

& — 


(13) 


m 


If the ving is flat, then for the same constant angle of attack over the 
wing for hoth direct and reverse flow 


v 0 ( S,n) = w 0 '(x,y) = -Ya 

over the entire wing. Tnus, for this caBe, the elementary lifts 8L 
and 8L' in equations (13) and (14), respectively, hecome identical. 
In hoth direct flow and reverse flow, the total wing lift is given by 
the integration of the elementary lift over the plan form, the element 
of integration being dx dy. Thus, the derivative for the wing 

in direct flow and reverse flow is 


°La 



The theorem of reversibility is thus proved for the lift of wings of the 
general type illustrated in figure 1 or 4. 

Derivatives Cm^ and Cl^.- If the wing in direct flow (fig. 4(a)) 
is at a constant angle of attack a, then 

v 0 (!-,tO = -Va 


(15) 



NACA TW 1943 


15 


The pitching moment 8M about the 7 -axis induced hy the sources in 
region A on the element dx dy is obtained by means of equation (13) 
and the value for v 0 given by equation ( 15) • Thus, 

8M(dx dy) = -x6L 



Wow if the wing in reverse flow (fig- 4(b)) is pitching about the y-axis 
with an angular velocity of q, then 

v 0 'U,y) = q* (17) 

(q is taken positive with respect to reverse-flow direction). The 
incremental lift 6L 1 on region A due to the sources in the ele- 
ment dx dy in reverse flow contributed by this pitching angular 
velocity is obtained from equation (l4) by substituting for w 0 ' the 
relation given by equation (17)- Thus, 



(18) 

A comparison of equations (l6) and (18) shows that the elementary 
pitching moment 6M' differs from the elementary lift 6L ' by a 
constant factor of proportionality, or 

6L ' = JLXM 

Va- 
in direct flow, the total pitching moment is given by the integration 
of the elementary pitching moment over th9 plan form, the element of 
integration being dx dy; similarly, in reverse flow, the total lift is 
given by the integration of the elementary lift over the plan form, th9 
element of integration being dx dy. This integration results in 


L 1 = JLM 
Va 
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When this relationship is converted to the nondimensional form 


there results 


CL a ’ 



C Lq' = 20^ 


(19) 


It should "be noted in connection with the relationship expressed in 
equation ( 19 ), that the reference axis for in direct flow and the 

reference axis for q in reverse flow must he fixed with respect to 
the wing configuration. For example, if C m for the win£ in direct 
flow is obtained with respect to an axis passing through the leading 
edge of the tip chord, , then when the flow is reversed, the leading and 
trailing edges become interchanged • The corresponding axis for the wing 
in reverse flow for pitching motion q then passes through the 
trailing edge of the tip. 


Damping in pitch Cm^.- If the wing in direct flow (fig* Ma)) is 
pitching about the y-axis, then 


v 0 U,t|) = -q .1 


( 20 ) 


likewise, in reverse flow (fig* 4(b)), 

v 0 'U,y) = I* (21) 

In figure 4(a) for the wing in direct flow, the elementary 
pitching moment 6M induced by the sources in 'region A on the ele- 
ment dx dy can be obtained from equation ( 13 ) by inserting the value 
for w Q given by equation (20) and multiplying by the moment arm x* 
Thus, 

8M(dx dy) = -x8L(dx dy) 


= -2pVqx dx dy 


-£(x - t)dt dn 

t>3 


1 / 1 - B tan y 
itB'Vl + B tan 7 




-i dS 


+ B 


( 22 ) 
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In figure ^(b) for the wing in reverse flow, the elementary pitching 
moment EM' induced on region A by the sources in the element dx dy 
is obtained from equation ( 14) by substituting for w Q 1 the relation 
given by equation (21) and by inserting the local moment arm I under 
integral Bigns. The moment arm for the last term is I = x. Thus, 


EM' = 2pVqx dx dy 


1 CP £(x - |)dl dn 1 f . 1 - B tan y \ P I d£ x 
* JJ A R3 nB Vl + B tan yj J Q R " B 


(23) 


A comparison of equations (22) and (23) shows that the elementary 
pitching moments EM and EM' are equal. For both the wing and its 
reverse, the total pitching moment is given by the integration of the 
elementary pitching moments 94 and EM' over the plan form, the 
element of integration being dx dy. Thus, the derivative for the 

wing and its reverse are the same: 



/ s EM dx dy 

a|v2sc 



The theorem of reversibility is thus proved for the damping in pitch. 

As indicated previously for the reversibility relationship for 

and Cp^, it should be noted that the reference axes for the moments and 

pitching angular velocities for the wing and its reverse must be fixed 
with respect to the wing configuration. 


Derivatives Cj^ and C^.- The derivative is a nondimen - 

sional form of the rate of change of lift with forward speed 



* f s( vc La) 

- ° (2I,) 
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A preceding section that treated the derivative Cp^ has shown 
that CLq^ = Cx^. The remaining parameters in equation ( 24 ) are also 
identical for the wing and its reverse ) therefore, 

CLu' = CLu 


The derivative C mu is a nondimens ional form of the rate of change 
of pitching moment with forward speed: 


n _ (&) 

- p _ ® | 


2 ^Sc 


w=Constant 


IfC^ 850 ^ Constant. 

gVSc 


- a & ( vt w) 


(25) 


A preceding section that treated the derivatives and 

has shown that 


0 ■ 0L3 

l ma = ~ 2 ^ 


where the axis of rotation is fixed with respect to the wing configura- 
tion. All the remaining parameters in equation (25) are identical for 
the wing and its reverse • The derivative ' for the reversed wing 
can therefore be obtained from the following relationship: 


c <'il?( vpL J (25a) 

where Cj^ is a function of V through its dependence on Mach number. 
Here again Cm u ' and Cp^ refer to an axis of rotation which is fixed 
with respect to the wing configuration. 


Derivatives C]^, Cp) a , and Cp)^.- 
for the derivatives Cp^, and Cp 


The ’ reversibility 
require special 


relationships 

consideration 
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for a general wing of the type illustrated in figure 1 or figure 4. If 
the wing has subsonic side edges, the portions that are raked out will 
develop a suction force giving thrust along the side edge- When the 
flow is reversed, these portions of the side edge are raked in and 
become, in effect, a subsonic trailing edge, and, therefore, they can- 
not produce any suction force if the Kutta-Joukowski condition iB 
applied. Consequently, the theorem of reversibility cannot be demon- 
strated for these cases. 

However, the resultant force on a thin flat wing with a supersonic 
leading edge and streamwise tip acts normal to the surface as there is 
no suction at the leading edge. Thus, the forces in the x-direction 
arise solely from skin friction. On the basis that the skin friction 
is unchanged when the flight direction is reversed 

°V - PDu 

If the skin friction is assumed to be independent of a and q, the 

derivatives Ct, and Ct* , with respect to body axes, will be zero for 
cl q 

both the wing and its reverse. 


Damping in roll Cjp.- If the wing in direct fiow, (fig- 4(a)) 


is rolling about the x-axis, then 


V 0 (t,T)) = -PT) 

(26) 

and, in reverse flow, (fig. 4(b)) 


^o'^y) = py 

(27) 


In figure 4(a) for direct flow, the elementary rolling moment BLj. 
induced by the sources in region A on the element dx dy can be 
obtained from equation (13) by inserting the value for w 0 given by 
equation (26) and multiplying by the moment arm y. Thus, 

6L r (dx dy) = -yBL(dx dy) 


= -2pVpy dx dy 


1 rr -q(x - |)d| dn 
* JJ A r3 


1 /l - B tan 7 \ r — q dl y 
«B \1 + B tan y JJ Q R + B 


( 28 ) 
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In figure 4(b) for reverse flow, the elementary rolling moment BLy ' 
induced on region A by the sources in the element dx dy is obtained 
from equation (l4) by substituting for wo ' the relation given by 
equation (27 ) , and inserting the local moment arm q under the integral 
signs. The moment arm for the last term is q = y. Thus, 


6Lj. ' = 2pVpy dx dy 


H q(x - 5)d£ dn 
f A r3 


1 f 1 - B tan r q d£ _ y 
«B\1 + B tan yJJ R B 


(29) 


A comparison of equations (28) and (29) shows that the elementary 
rolling moments BLj. and BLr ' are equal. For both the wing and its 
reverse, the total rolling moment will be given .by integration of the 
elementary rolling moments BLr and 81^. ' over the plan form, the 

element of integration being dx dy. Thus, the derivatives for 

P 

the wing and its reverse are the same: 




51^ dx dy 




UNSTEADY MOTION ' 

Derivation of' Basic Equations for Unsteady Motion 

The velocity potential for unsteady motion at the point (x,y) in 
tho plane of the wing and at time t may he expressed in the following 
form (see reference 9 or reference 10): 


</(x,y,t) 


1 /T CwoU^t - ti) + w Q U,q,t - T2)]dg dq 

2jt JJ R 


whore T p and T 2 are the time-lag functions 


(30) 



(x - Dm _r_ 

2 2 
B a B a 


(3D 
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T o = (* . -J ) M - JL_ (32) 

2 B 2 a B 2 a 

and -where the integration extends over the entire region in the plane of 
the wing that can influence the point (x,y) . 

Equation (30), which gives the potential distribution with time 
dependency, corresponds in form to equation (l), which gives the 
potential distribution for the steady case. The significance of the 
form of equation ( 30 ) has been discussed previously following the 
presentation of equation ( l) . The region that can influence a 
point (x,y) in the unsteady case is contained within the Mach forecone” 
from the point, and corresponds identically to that described previously 
for the steady case, as indicated in figure 2. 

The present analysis considers unsteady motions in which the wing 
is accelerating uniformly in the vertical direction. It is shown in 
reference 9 that, for this type of unsteady motion in certain cases, the 
influence of the external field F (fig. 2) on the point (x,y) is 
identically canceled by the effect of the region G in the wing. This 
cancellation is valid under the conditions noted previously for the 
steady case. 

The angle of attack at the beginning of the unsteady motion is 
assumed to be zero, inasmuch as this assumption does not result in any 
loss in generality because the principle of superposition is applicable. 
Then, for a uniform acceleration, w 0 (£,ii,t) = -Vat for all values of | 
and tj. Thus, 


■WoU^t - ti) = -Va(t - ti) 


v 0 U,Tj,t " t 2 ) = -Va(t - t 2 ) 


If these values for w 0 are substituted into equation (30), the 
following equation results: 


>t) = I 


/ 

in 

a 

, M(x__- |1 

V 

t' B^q j 

11 UK 

j 



d£ dtj 


( 33 ) 


The pressure difference at the point (x,y) at time t is 
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The quantity (jf x may be obtained by differentiating equation (33) 
under the integral sign with special attention to the variable l imi ts 
for I and q. (See fig. 2.) If this differentiation is performed, 
the contribution of each component region for ^(x,y,t) in equation ( 33 ) 
results essentially in five integrals which are analogous to the result 
obtained for the steady case. (See also appendix A for steady case.) 
The quantity 9^t niay be obtained by direct differentiation of equa- 
tion ( 33 ) under the integral sign. 


As a result of an analysis similar to that indicated in appendix A, 
the pressure difference from equations ( 33 ) and (34) can be expressed 
in the following form: 


/ 


£P(x,y,t) = 2pV 2 ct 


V 


pp 

t M(x - |) 

(x - £)d! dq 

t - ‘"vy 

L B 2 a J 

JJ A R 3 



n 


JU A 


dl dn 
R 



■where A and a r,efer to the area and line, respectively, as shown in 
figure 2 . Letting t— >0 in equation ( 35 ) in order to obtain the 
instantaneous value of the derivative and simplifying yields 


AP(x,y) 


2pV& 

M 2 

rr 

(x - |) 2 d| dn 

«B 2 


A 

R 3 



(x - £)d£ IT dl dn 
R ‘ JJ A R 


(36) 


Elementary Lifts for Wing in Direct Flow and in Reverse Flow 

Equation ( 36 ), when modified for the wake effect contributed by 
portions of a subsonic side edge that are raked inward, is the basic 
equation for deriving a relation between the lift contributed by a 
uniform variation of angle of attack with time for the wing and its 
reverse. The concept of a forezone of influence and an afterzone of 
influence for a point is, again used in the same manner as described 
previously for figure 3 for the case of steady motions. 
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. The following results are obtained by means of a development 
analogous to that used to obtain the elementary lifts for steady motions 
given previously. 

Elementary lift in direct flow .- The lift on the element dx dy in 
figure 4(a) for direct flow can be evaluated at once. The pressure 
difference at dx dy is given by equation ( 36 ) and the lift 5L is 
equal to AP dx dy. Thus, since 

= . l/ l - B tan y \ 
dx B\1 + B tan 7 / 


&L(dx dy) = dx dy 

«B 2 


M c 


H 


(x - I ) 2 dS dq 
B3 


+ $£ ( 1 - B tan £ ] r (X - S.)dl . rr dl dn 
B \1 + B tan yj J a B JJ A B 


(37) 


Elementary lift in reverse flow .- The lift on region A of fig- 
ure 4(b) due to the sources in the element dx dy in reverse flow is 
less easily obtained. As noted previously in the analysis for steady 
motions in connection with the reverse flow shown in figure 4 (b), 
account must be taken of the vorticity shed in the wake of the portion 
of the wing tip that acts as a subsonic trailing -edge. (See fig. 3(b).) 
The contribution of this vorticity or wake effect to AP ’ for a point 
on the wing in reverse flow is designated as A(AP'). This contribu- 
tion A(AP’) is determined on the basis that the Kutta-Joukowski condi- 
tion is satisfied along the subsonic side edge. The contribution 
to AP’ of the part other than the wake effect is obtained from 
equation ( 36 ) with reference to figure 3 (a) in a manner analogous to 
that used in the steady -motion case to derive equation ( 4 ) from 
equation ( 3 ). By following this procedure, AP ’ at the point (6,tj) in 
reverse flow (fig. 3 (a)) is obtained as 


ap’u,ti) = zm 


jtB 




J} 


T (x - g) 2 dx dy 
A r3 


/ 1 + B tan 


B \ 1 - B tan' 


il) [ i*,.- $)ax . IT a* dr 

: 7/ J a' B JJa B 


+ A(AP ') ( 38 ) 
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An analysis for determining the wake -effect correction for uniformly 
accelerated motion A(AP') is presented in appendix C. As a result of 
this analysis, it can be shovn that, analogous to the case of steady 
motions for flat plates in which the side edges are subsonic and raked 
in, if the Kutta-Joukowski condition is applied, the vorticity in the 
wake in uniformly accelerated motion is of just sufficient strength to 
cancel the contribution to AP of a line integral vhich yields an 
infinite pressure along the subsonic edge. This line integral which is 
canceled by the wake effect is evaluated on the wing along the line a ' . 
As a result of the analysis in appendix C, the wake -effect contribution 
to the pressure difference in equation ( 38 ) is found to be 


A(AP') = - ^P v< * m 2 tan L r (* - lid* 
«B 2 1 - B tan 7J a ' R 


(39) 


If the expression for A(AP') in equation (39) is substituted into 
equation ( 38 ) and the result is simplified, the expression for AP ' at 
point (!,tj) in figure 3 (a) in reverse flow becomes 


£P'(£,t]) 


2pVd 

jtB 2 



(x - g) 2 dx dy 
R3 


+ * r (* - ildx _ rr dx dy 

B J a 1 B JJ A R 


m 


The lift on region A due to the sources in the element dx dy 
in reverse flow (fig. 4(b)) is obtained from equation (40) in a manner 
analogous to that used in the analysis for steady motions in deriving 
equation (12) from equation ( 6 ) . Thus, the differential components of 
equation (40) are interpreted as elementary lifts induced by particular 
source elements, and these lifts are then resummed in a different 
fashion to give the lift induced on the region A by the sources 
in dx dy. The contribution of the two surface integrals in equa- 
tion (40) can be obtained in a manner similar to that indicated in the 
preceding analysis for steady motions, and the contribution of the line 
integral can be obtained in a manner analogous to that described in 
appendix B. This procedure leads to the result that the total lift 
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induced on the region A by the sources in the element dx dy in 
reverse flow is given by the following expression: 



Reversibility Relation for Lift in Unsteady Motion 

A comparison of equations (37) and (4l) indicates that they are 
identical; in both direct flow and reverse flow, the total wing lift' 
will be given by the integration of the elementary lift over the plan 
form, the element of integration being dx dy. Thus, the deriva- 
tive for the wing in direct flow and in reverse flow is 


°Ld = 


The theorem of reversibility is thus proved for the lift contributed by 
a uniform variation of angle of attack with time for wings of the 
general type illustrated in figure 1 or figure 4. 


RESULTS 


The reversibility relations that have been derived in this analysis 
were obtained for a subsonic wing-tip configuration of the type shown 
in figure 1 or figure 4. All the results obtained for this illustrative 
plan form will likewise apply to wings with streamwLse tips or those 
with subsonic tips of arbitrary rake as long as the edges are nonre- 
entrant. Furthermore, as indicated in the preceding analysis, equa- 
tion ( 2 ) for steady motion and equation ( 36 ) for unsteady motion may 
be applied also to any wing having supersonic tips of arbitrary rake 
since the component regions C, D, and E simply degenerate into 
regions of the type C and E. Therefore, the relations obtained 
between a wing and its reverse for the stability derivatives for the 
case of the streamwise tip will apply al3o to wings that have super- 
sonic tips of arbitrary rake. 
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A complete summary of the relations between the characteristics of 
a wing and its reverse in supersonic flow is presented in table I. The 
table includes all the results obtained in the present paper and also 
those obtained by previous investigators. It should be noted that all 
the stability derivatives refer to a system of body axes which is fixed 
with respect to the wing configuration. 


CONCLUSIONS 


An investigation has been made by means of the linearized theory 
of the relations between the supersonic stability derivatives of a wing 
and its reverse obtained by interchanging the leading and trailing 
edges. Th 9 analysis included wings that have a supersonic leading 
edge, a supersonic trailing edge, and either supersonic side edges or 
nonre-entrant subsonic side edges that are noninteracting, that is, . 
neither edge lies in the zone of action of the other. The investigation 
was made for thin flat -plate wings -and included steady vertical apd 
longitudinal motions, steady rolling and pitching, and motions in 
which the wing is accelerating uniformly in the vertical direction. 

The results of the investigation were obtained with respect to body 
axes and indicated the following significant conclusions: 

1. The rates of change of lift and drag with angle of attack and 
forward speed, the damping in roll, and the rate of change of lift with 
a uniform variation in angle of attack with time were all equal for the 
wing and its reverse . 

2. With respect to axes of rotation that were fixed with respect 
to the wing configuration, the damping in pitch was equal for the wing 
and its reverse j whereas, the rate of change of pitching moment with 
angle of attack for the wing was equal to the rate of change of lift 
with pitching velocity multiplied by the stream velocity for the reverse 

of the wing = Jd- 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., June 15, 19^9 
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APPENDIX A 

DERIVATION OF EXPRESSION FOR DISTRIBUTION 


The linearized velocity potential for steady motion is from 
equation (l) 


[x,y) = - i 

It 


W 0 (l,Tl)^ dl l 


]f(x - |) 2 - B 2 (y - tj) 2 


(Al) 


where the integration is to be performed over the entire area that 
contains singularities that can influence the point (x,y) . The regions 
of influence for a point (x,y) are contained within the Mach forecone 
from the point and are illustrated in figure 2. In figure 2, the wing 
leading edge is of the supersonic type) whereas the wing tip is of the 
subsonic type* The regions C, D, E, F, and G contain singularities 
that will influence the point ( x, y ) . In the case of region F, these 
singularities are located off the wing tip. It will be noted in 
figure 2 that, in general, the limits of integration for • I and T| for 
a particular region may be functions of the variables x and y. 


The quantity ^ x (x,y) may be obtained by differentiation of equa- 
tion (Al) under the integral sign, with special attention to the 
variable limits . Thus, for an arbitrary region of integration, the 
following expression is obtained: 

12 

w 0 ( |,T|)(x - |)d | dq 

W 


?xU,y) = ' 



[U - t)* - Is (y 
v Q ( ^ T l 1 )dl 


x - |) 2 - B 2 (y - q x ) 2 


^p(l 

/(x - 5) 2 - b 2 (j - ne) 2 


>>3 


(continued on next page) 
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" ST" 


p ^2^*^ * ^1^ 


U 


tJi(x,y, ij) 


y 0 ( 6 i»n)^ 

V(x - I x ) 2 - B 2 (y - ri ) 2 


dl Q PVr(*>*>Z2) 

+ 


Wo(^2,Tl)<iTl 


dx 


Tii(^>y> 62) 


|/(x - i 2 ) 2 " B 2 (y - *i)‘ 


> 


(A2) 


In the applications of equation (A2), the integrals become improper 
■when the region of integration is bounded by any part of the Mach fore- 
cone from the point under consideration. In such cases the procedure 
is to alter the l imi ts of integration for t) and | by means of a 
limiting process. Thus, the limits t^, £]_, or Zq become rj^ + e, 

T 12 _€ >^l + € > or ^2 _€ * The final result is then obtained by 
permitting 6 to approach zero. 

Equation (A2), as expressed for each region of integration, con- 
sists of one surface integral and four line integrals evaluated around 
the boundary of the region of integration. Consider the regions- C, D, 
E, F, and G in figure 2. It has been indicated in the body of the 
paper that; for a thin flat plate, the effect of region F will be 
identically canceled by the effect of region G. (See discussion 
following equation (l)«). Therefore, equation (A2) can be considered 
only with respect to the regions C, D, and E. 

/ 

The subsequent analysis refers to equation (A2) and to figure 2. 
This equation and this figure have been set up for the condition that 
the integrations are always performed first with respect to q and 
then with respect to |. Lines are indicated by means of subscripts 
that refer to the limit and to the region of integration; thus, 

indicates the line along the upper limit for the C region of 
integration (fig. 2). 
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It is apparent from equation (A2) and figure 2 that the following 
line integrals are identical and will cancel each other: 


d^2C . 
“§T- 


d§2D 

^x 


dTj = - 


hi 


ID 


dx 


B 2C 


=2D 


• • • 


, J-Ci 

a i * - sr 


Furthermore, the integral 


’ID 


<^IE 

,X 'llE 


dn 


dt) 


> 


y 


(A3) 


dx 


. dr) = 0 


’1C 


■because the upper and lower limits are identical. The integral 


l 2E = 


ai 2E p l te( x >y>*2E) 


v o(^2E^J dT l 
l/( x - 5 2E^ - ® 2 (y - D‘ 


must be evaluated as the limiting value of the integral as 

^2E * x 


(A4) 


and 


This integration yields 


►y. 


I 2E^2E 




v 0 (x,y)it 
y = — 


Inasmuch as tiqc 



independent of 




dx 



dl = 0 


(A5) 


(A6) 
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The integrals (see equation A2) 
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J 


(A7) 


are evaluated along boundaries of the Mach forecone from the point (x,y) . 
Along these boundaries, the denominators of the integrands of these 
expressions are zero. Thus, each of the integrals in the expres- 
sions (A7) gives infinity. However, the surface integral in equa- 
tion (A2) also has as limits 


^lC' ’liD' ^lE' an & T l2E 


(AB) 


If the surface integral in equation (A2) is evaluated for regions C, D, 
and E for conditions in which v 0 is either constant, is any function 
of |, or varies linearly with tj, the substitution of the limits given 
by (A8) will give infinities which identically cancel the infinite line 
integrals given by the expressions (A7)« 

The foregoing analysis thus shows that the evaluation of' $ x (x,y) 
is obtained from equation (A2) by Ineans of the finite part of the 
surface integral evaluated over the regions C, D, and E, or A by 
the expression given by equation (A5) and by the line integral evaluated, 
over Tjgp or a. Thus, 


0 x ( x >y) = Z 


v 0 ( 1,0 (* - Qd-I dT ) 

R 3 


i ( P v 0 (i,n)di w o( x ,y ) 

■ *\-£7j a s B 


(A9) 


A further discussion of equation (A9) is given in the body of the 
paper immediately following the presentation of equation (2). 
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APPENDIX B 


EVALUATION OF CONTRIBUTION TO LIFT ON REGION A OF LINE 
INTEGRAL ALONG a ' IN EQUATION (6) IN REVERSE FLOW 


The contribution to the lift on region A (fig- 4(b)) of the line 

1 /’ w n '(x, 

integral — / — 

itB.j a , R 

special treatment. It is convenient to transform the integral to 
.the u,v Mach line coordinate system shown in figure 6. In this 
coordinate system 


y)dx 


■which appears in equation (6) requires 


A 


» - ^ - B i); v = s (l + B,l) 


! 1 * *)i 1 = r(' ■ "I 

M M 


(Bl) 


The element of area in the u,v coordinate system is given by 


j/liH'',du dv = du dv 
Vu,v/ m2 

The u,v axes coincide with the Mach lines from the origin. The 
line a' is a Mach line (see fig- 3(a)) that can be specified by the 
relation 


v = a 1 = Constant 


(B2) 


(More rigorously, the equation for a' should specify limitations 
for I, but this specification is not needed in what follows.) In the 
reverse flow indicated in figure 6, the field element is denoted 
by du dv and the source element by dvq du^ . Along the line vj_ = a 

dv^ = 0 and, by equations (Bl), dx = ® du^« Thus, the line integral 


1_ P v o'Ulj a 

*MJa' R 

is to be regarded as expressed in terms of u, U]_, v. 


(B3) 


where R 
and a ' . 
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Refer now to figures 3(a) and 4(b) for reverse flow. In this 
reverse flow, let the line a in figure 4(b) pass through a point 
corresponding to P shown in figure 3(a)* Then the forezone from a 
point lying along the line a in figure 4(b) will b9 identical with 
the forezone A' from the corresponding point P in figure 3(a). 

Thus, in figure 4(b), the edge v of the afterzone from a point (x,y) 
lying on v, in the source element dx dy, corresponds to the edge a' 
of the same after zone which may be drawn from a point on a ' in 
figure 3(a). The point (x,y) in the source element dx dy of 
figure 4(b) then may be considered as lying along a' in figure 3 (a). 
This situation is represented in Mach coordinates in figure 6. In Mach 

OD 

coordinates, x,y becomes uq,vq. Increments of area — -2. du dv 

M 2 

pT> 

and — du]_ dv^ are shown at u,v and u^, V]_, respectively. Equa- 
M^ 

tion (B2) specifies line a' and line a is specified by 

u = a (B4) 


On the basis of the foregoing considerations, the line integral (B3) 
.gives the velocity Uq ' at point P contributed ultimately by all the 
sources along the line a'. Those sources in the length duq will 
contribute the part 


62Uo'(u,v) 


_1_ w o'( u l> a, ) du l 
R 


(B5) 


Note that the contribution at an interval dv removed from point P 
along the line -a differs from that given by equation (B5) by an 
infinitesimal of higher order. Now consider the contribution at an 
interval du removed from P along the other Mach line from P 
(parallel to a'); the sources along a' do not contribute, but those 
along a line parallel to a ' and displaced from it by an internal dv 
do contribute. Again the contribution differs from that given by equa- 
tion (B5) by an infinitesimal of higher order. Accordingly, equa- 
tion (B5) may be taken to give the contribution to the velocity Uq ' 
anywhere within the element du dv due to the sources in the 
element duq dvq. 

The increment from the line integral to the lift on the 
element du dv is 2pV^^)du dv times the velocity Uq ' given by 
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equation (B5) • The incremental lift on all the elements du dv 
touching the line a (fig. 6) is thu 3 


V- ■ sfgy-,.- 

In the present case, a' is simply vp, and, as can be obtained from 
figure 6 

fl - B tan 

4u = vrmsr7/ d 1 


Thus 


r^t > _ 2pV/2B\/l - B tan y\, , 

^ ' mitfJ KTTT uZ 7) ^1 


w o '( u l> v l) 4v 
B 


If this equation is converted to Cartesian coordinates by means of 

equations (Bl) , and since — dv = d£ along a, there results 

M 


82L' 


2pv/l - B tan y\ 
«B \1 + B tan 7) 


dx dy 


v 0 '(x,y)di 

R 


(B6) 


Equation (B6) expresses the second part of the incremental lift contri- 
buted to the region A of figure 4(b) by the sources in dx dy. This 
second part results from the line-integral term in equation (6) and, as 
has been seen, is localized along the line a. 
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EVALUATION OF WAKE -EFFECT CONTRIBUTION TO PRESSURE BASED 
ON KDTTA-JOUKDWSKI CONDITION ALONG SUBSONIC RAKED -IN 
SIDE EDGES IN UNIFORMLY ACCELERATED VERTICAL MOTION 


A solution for the contribution of the wake effect that satisfies 
the Kutta-Joukowaki condition or A(AP') in equation ( 38 ) must meet 
the following requirements: 

1. The resultant pressure difference AP 1 in equation ( 38 ) must 
vanish along the subsonic raked- in portions of the side edge and at all 
points exterior to the wing. 

2. The solution A(AP') must give zero vertical velocity at all 
points on the wing surface . 

3* The additional potential associated with A(AP') must 
satisfy the linearized partial-differential equation for unsteady flow, 
that is. 


An: + 


yy 


+ t 


ZZ 


2M 2 * 
— y xt 


- = 
«2 ¥ tt 


(Cl) 


In the analysis to follow, it will be convenient to deal at first 
with a wing in direct flow as shown in figure 4(a), but the side edge 
will be assumed to be raked in instead of raked out. Thus, the field 
point will be denoted by (x,y) and the source points by (£,tj). 

A consideration of the pressure difference obtained frcm 
equations ( 35 ) and ( 36 ) shows that AP will be made to vanish along 
the subsonic side edge and at all points exterior to the wing by the 
addition of the term 


A(AP) 


4pV^d tan 7 
it(l + B tan 7) 



5) 


dl 


(C2) 


In equation (C2), the integral is to be evaluated along a line 
corresponding to a in figure 4(a), and 7 is the local slope of the 
side edge at the point of intersection of line a with the side 
edge (-7 in fig. 4(a)). The solution given by equation (C2) fulfills 
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the first requirement noted previously for the vanishing of the 
resultant pressure difference along the subsonic edge and off the wing. 
It remains, therefore, to show that requirements 2 and 3 are fulfilled 
by the solution (C2) . 

Equation (C2) may be expressed in the form 


A(AP) = K 




M^x \ r dl 

^Y/Ja R 


(C3) 


where 

K = ^PV^d tan 7 
Jt(l + B tan 7 ) 

Consider the first integral in equation (C 3 ) . This integral is 
identical -with the integral -which gives the -wake -effect contribution 
for a steadily pitching wing -when the Kutta-Joukovski condition is 
applied along the subsonic raked-in portions of the side edge, as 
expressed by equation ( 5 ). It vill be noted that the notation in 
equation ( 5 ) is for reverse flow where the field point is denoted 
by ( I, tj) and the source points, by (x,y). Thus, for a wing pitching 
steadily at unit velocity, w 0 ' = x in equation ( 5 ) corresponds 
to w 0 = I in the first integral in equation (C3) • In order for 
equation (5) to represent a valid solution for the wake effect, the 
potential associated with this Auq ' contribution must give zero 
vertical velocity on the wing surface and must satisfy the partial- 
differential equation for steady motion. Therefore, the potential 
associated with the first integral in equation (C3) likewise gives 
zero vertical velocity on the wing surface and is a solution for the 
partial-differential equation, for steady flow. If this potential 
satisfies the partial-differential equation for steady flow, 

p 

then -B + tyyy + ijr zz = 0} also, since this flow is independent of 

time, the terms ih equation (Cl) that involve and are zero. 

Thus, the potential i|r satisfies the unsteady partial-differential 
equation (Cl) . 

Consider now the second term in equation (C3) • The 

integral J* ^ is identical with the integral that gives the wake- 

effect contribution for a wing at a constant angle of attack when the 
Kutta -Joukowski condition is applied along the subsonic raked-in 
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portions of the side edge, as expressed by equation (5) • For this 
case, equation ( 5 ) assumes the form 

Au » F 7 M (C4) 

J a R 

where F is a constant. In order that the Au in equation (C4) be a 
valid solution for the vake effect, the potential associated with it 
must satisfy the partial-differential equation for steady flow and must 
give zero vertical velocity on the wing surface- 

The foregoing conclusions regarding equation (C4) will be used., to 
show that the potential associated with the pressure contribution of 
the second integral in equation (C3) will satisfy the partial- 
differential equation for unsteady motion (Cl) and will also give zero 
vertical velocity on the wing surface. Let this potential be denoted 
by \|r and let the potential associated with the pressure contribution 
given by equation (C4) be denoted by t. The AP contribution of the 
second integral in equation (C3) is 

. (C5) 

Since t x = F / — (see equation (C4)) 

• J a F 



Substituting this relationship into equation (C5) results in the 
equation 



(C6) 


K 

~ 2pVF 


where 
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Integration of equation (C6) yields 

(* [* n <°7) 

*11!^ ’ ^ X LE 

(■where x 1 is an au xi l i ary variable of integration which replaces x) 
and integration of equation (C7) by parts results in 




(C8) 


Differentiation of equation 


(C8) yields the following results: 


*xx «' E 


t - 


M^x 


W 


~2~ ]' xx ~o~ T X 

B^vy B^V 


V^tx 


*: = E 

T yy 


i | r = E 

zz 


t - 


M c x 


-V^r + 


w 


B^V/ B^Y 


T yy 4x1 


2 le 


I f *yyt 

^ *LE 


dx' 


> (C9) 




-B C V 


t dx' 

ZZ 


C LE 


1 

V 


*zzt ta ' 


^LE 


*xt = ET X - br- 


tt 


In the present analysis only a linear variation of ' a with time 
is considered; consequently, the potentials associated with this linear 
variation of a. will also vary linearly with t. (See equation (33).) 
Thus, in this analysis, the terms that involve the second derivative 
of with respect to t ■ will be discarded. The partial-differential 
equation for unsteady motion, equation (Cl), then becomes 


-B^xx + 


+ 

yy + 


zz 



/ (CIO) 
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Substituting the relations given by equation (C9) into the left-hand 
side of equation (CIO) and rearranging the terms results in 



Expression (Cll) may be simplified by use of the fact that, as indicated 
previously, t satisfies the partial-differential equation for steady 
motion so that 


-® 2t xx + T yy + T zz = °. (C12) 

.If the relation given by equation (CI2) is substituted into the 
expression (Cll), all the terms in expression (Cll) will vanish with 
the exception of the term 


1 

y 


r ( 

'-' X LE V 


1 fcnct + tyyt + ^zzt/^' 


It is seen from equation (CIO) that 


-B 2 ^ + * + t _ 2M 2 - 

xx + "yy + \z ~ — ^xt 


therefore 


^xxb + ^yyt + ^zzNt = “v“^xtt = ® 


(C13) 
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Thus, the expression (Cll) vanishes, and the potential \|r satisfies 
the partial-differential equation for unsteady motion, equation (Cl) . 

The contribution of the potential ijr to the vertical velocity is 
obtained by differentiating equation (C8) with respect to z. Thus, 


t 


0 






1 

V 



dx' 


(Cl4) 


Differentiating equation (,C8) first with respect to t and then with 
respect to z yields 


*t = Et 

and 


*tz =■ ET z 


(C15) 


It has been shown previously that the vertical velocity contributed 
by r on the wing surface is zero, or 



(C16) 


Then, substituting the relations (C15) and (Cl6) into equation (Cl4) 
results in 


t,_,o " 0 ■ ' (C17) 

The expression given by equation (C2) therefore satisfies all the 
requirements for a valid solution. Conversion of this expression so 
as to apply to the reverse flow shown in figure 4(b) yields 

A( AP ' ) = - 4pVaM 2 tan 7 

itB^(l - B tan 7) 

which is equation (39) given in the body of the paper- 


(x - |)dx 
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TABLE I 

SUMMARY OF RESULTS FOR RELATIONS BETWEEN CHARACTERISTICS OF A WING AND 


ITS REVERSE IN SUPERSONIC FLOW 


[a 11 results are vith respect to body axes fixed in wing; results for 
subsonic tips are based on application of Kutta-Joukovski 
condition^] 


Characteristics or 
stability derivative 

Relation between 
wing and its reverse 

Plan -form 
configurations 
(c) 

Source 

°»vo 

Ctj * = Cjj 
^ vo , v ^wo 

(a) 

All plan forms 

References 1 to 4 


V " Cl * 

A 

References 2 to 4 
and present paper 

B 

Reference 4 and 
present paper 

C “a 

v-%* 

0>) 

A 

Reference 3 

B 

Present paper 

c Da 

11 

? 

II 

O 

A, B X 

Present paper 

CL a 

<V = ^ 
0>) 

A 

Reference 3 

B 

Present paper 

C niq 

c m^' = c m< 1 ‘ 
0>) 

A 

Reference 3 

B 

Present paper 


CD 4 ' = CD q 
(b) 

A, Bj_ 

Present paper 

°Lu 

<V = Cl„ 

A, B 

Present paper 

Cmu 

(b) 

A, B 

Present paper 

^Du 

QDu' = °I>U 

A, B X 

Present paper 

Clp 

^Zp 1 = CZp 

A 

Reference 3 

B 

Present paper 


c l*’ = 

A, B 

Present paper 


a If leading edge is subsonic, airfoil section must be symmetrical. 


^Axes for rotations for wing and its reverse are fixed vith respect to ving configuration. 

°The notations A and B have the following significance: 

A all edges are supersonic; supersonic side edges are raked arbitrarily 
B supersonic leading edge, supersonic trailing edge, subsonic side edges, are nonre- 
entrant and noninteracting, that is, neither edge lies in the zone of action of 
the other; when used- with the subscript 1, plan form is restricted to Btreamvise 
tips (see fig. 1). 
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Figure 1.— General illustrative plan form for analysis. Subsonic edges do not interact. 
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-Leading edge^ 
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. » T) 

(x y y ) or P 


Regions C+D+E=A- 
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/ 

/ E 


y 

'E^Mach angle 


■Vlmgtip 


/ (xy)orQ 


Limits of Integration When Integrations are reformed 
Erst with Respect to no and Then with Respect to g 



Figure 2.— Regions of influence for a point on a wing. Supersonic 
leading and trailing edges with streamwise wing tip. 













Figure Positive directions for moments and angular velocities in direct flow and in reverse flow. 
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x or g 


Figure 6.- Mach line coordinate system. 
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